Transports of Disordered Carbon Nanotubes with Long Range Coulomb Interaction 
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Transport properties of disordered carbon nanotubes are investigated with including long range 
Coulomb interactions. The resistivity and optical conductivity are calculated by using the memory 
functional method. In addition, the effect of localization is taken into account by use of the renor- 
malization group analysis and it is shown that the backward scattering of the intra-valley and that 
of the inter-valley cannot coexist in the localized regime. 
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Single wall carbon nanotubes (SWNTs) are the new 
materials which are an experimental realization of one- 
dimensional (1-D) electron systems [Q. Since the SWNT 
is made by rolling up a graphite sheet, it is expected that 
fascinating properties different from the conventional 
quantum wires made from semiconductor heterostruc- 
tures will be observed. From this point of view, transport 
properties of disordered SWNTs have been discussed in 
Refs. [|| and ||. It has been predicted that the back- 
ward scattering due to the impurities vanishes when the 
range of the impurity potential is much larger than the 
lattice constant, but that the backward scattering reap- 
pears when applying a magnetic field perpendicular to 
the tube axis. It should be noted that the absence of the 
backward scattering holds for the graphite sheet as well 
as any SWNTs. 

In the above studies, electron correlations have been 
neglected. The 1-D nature together with the electron- 
electron interaction has been known to result in a variety 
of correlation effects in SWNTs in case of short range in- 
teractions 0,5],[§ and for the long range Coulomb interac- 
tion [0,||||[lO|. Effects of electronic correlation in SWNTs 
has been measured in the Coulomb blockade regime |Tl| l 
as well as for Ohmic contacts [^2j. In the latter exper- 
iment, power-law dependences of the conductance as a 
function of temperature and of the differential conduc- 
tance as a function of bias voltage have been observed 
and interpreted in terms of tunneling into clean SWNT 
with the interaction. 

Even when the Coulomb interaction is taken into ac- 
count, the conclusion of the absence of the backward scat- 
tering in Refs. |2) and j| is not changed. However, the 
effects of the Coulomb interaction on the transport in dis- 
ordered SWNTs should be observable in case of shorter 
range impurity potentials. In the present paper, I will 
discuss the transport properties of SWNTs with short 
range impurity potentials and the long range Coulomb 
interaction. It is shown that the interaction gives rise to 
an enhanced resistivity compared to that without the in- 
teraction. In addition, the interaction leads to a power- 
law dependence of the resistivity as a function of tem- 
perature and modifies the power of the frequency for the 
optical conductivity. In the localized regime, I find that 



intra-valley and inter-valley backward scattering cannot 
coexists. 

The SWNT has metallic bands when the wrapping 
vector, w = N + a + + iV_a_, satisfies the condition, 
N + - N_ =0 mod 3, where a± = (±l, v / 3)a/2 with 
a being the lattice constant. The Hamiltonian of the 
metallic SWNT with the long range Coulomb interac- 
tion is written by the slowly varying Fermi field, ip pa s, 
of the sublattice p = ±, the spin s = ± and the valley 
a = ±, as follows (l(J, 

Ho = -ivo ^2 ae~ ipax J dxipl a8 d x ip- pas 



V(0) 



dxp(x) 2 , 



(1) 



where p(x) = Y^pas^pas^Ppas, "o - 8 x 10 5 ?7i/s, and 
V(0) = 2e 2 /Kln{R s /R) with k ~ 1.4, R s and R be- 
ing the cut-off of long range of the Coulomb interac- 
tion and the radius of the tube, respectively. In Eq. ([!]), 
X = tan -1 [(iV + - iV_)/V3(iV + + N-)] is the chiral angle 
(x = corresponds to an armchair nanotube). In the 
above expression, I disregard the matrix elements of the 
interaction of the order of a/R, which lead to energy gaps 
0,||J|Jl^]. Therefore, the present theory is valid when the 
temperature, T, or the frequency, w, are larger than the 
gaps induced by the Coulomb interaction. 

The impurity potential introduced as disorder of the 
atomic potential is given by the following Hamilto- 
nian (|, U lmp = ^paa's I dxV L'( x )^ pa s^ pa 's, where 
V? ; (x) is the impurity potential at the sublattice, 
p, by which the electron on the valley, a', is scat- 
tered into the valley, a [^3|. Here I diagonalize 
the kinetic term in Eq.(nl) and move to the basis of 
the right-moving (r = +) and the left-moving (r = 
— ) electrons by the unitary transformation, ip pas = 



Then, Hi mp is writ- 



( e -ipax/2/^2) ^ r ( ar )( 1 -f)/ 2 V'r 
ten as follows, 

Himp = 2 / dX { i V ^( x ) + V aa( x )] ^tasAas 
ras ^ 

+ W ax V+- a (x) - e- ia W-_ a (x)] CsVV-a, 
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+ [V2a(x) - V- a (x)] tUji-ra. 

+ [e iax V+_ a (x)+e- ia W-_ a (x)} ^ Q8 V- 



(2) 



Here the first (second) term expresses the intra-valley 
(inter-valley) forward scattering, and the third (fourth) 
one is the intra-valley (inter-valley) backward scatter- 
ing. When the range of the impurity potential is much 
larger than the lattice constant, V+ a — V~ a and V^_ a = 
V~_ a — 0, which leads to vanishing of the third and 
fourth terms in Eq.(||), i.e., the absence of backward 
scattering, in agreement with Ref. 11). Here I consider 
the case of an impurity potential with range shorter than 
the lattice constant by retaining finite matrix elements 
of the backward scattering in Eq. (^J) . I disregard forward 
scattering because it does not contribute to transport. 

As was pointed out by Abrikosov and Ryzhkin Jli) , 
in 1-D systems and in the limit of weak impurity 
potentials, the interaction between the electrons and 
the impurities can be parameterized by uncorrelated 
Gaussian random fields. Here I extend this method 
to the present model and introduce two kinds of the 
random fields, 7](x) and £(x), expressing the intra- 
valley and the inter-valley backward scattering, re- 
spectively. The fields have the probability distribu- 
tions, P v = exp { — (2Z>i) _1 J dx-ri 2 (x)} and = 
exp { — (_D 2 ) -1 f dxt;(x)£,* (x)} where D x and D 2 are 
given by vq/ti and vq/t2 with n (r 2 ) being the scatter- 
ing time due to the intra-valley (inter-valley) backward 
scattering. The Hamiltonian of the impurity potential is 
given by 



h.c. 



(3) 



Note that the intra-valley (inter-valley) backward scat- 
tering, where the momentum transfer in the scattering 
process is small (large), is parameterized by a real (com- 
plex) field. 

Here I utilize the bosonization method and intro- 
duce the phase variables expressing the symmetric (8 — 
+) and antisymmetric (8 = — ) modes of the charge 
(j = p) and spin (j = a) excitations, 9jg and 
4>j$. The phase fields satisfy the commutation relation, 



[0js(x), (f>j>s>(x')} = i(7r/2)sign(a; - x')8j 
field, ip ras , is expressed as 



The Fermi 



where 8 n 



r 

: exp irq F x + — {6 as + r(p as \ 



(4) 



s6„> + a0„ 



Jas — "p+ r o\7(j -|- ~r~ mu p — olsO (j — and (j> as — (j> p + + 
■s(f>a+ + ct<fip- + as(j) a - . Klein factors rj ras are introduced 
to ensure correct anticommutation rules for different 
species r,a,s, and satisfy [rjras, r]r'a>s>}+ = 28 rr '8 a a'8 SS ' . 



The spin-conserving products r] ras r] r i a i s in the Hamilto- 
nian can be represented as Jt|, A ++ (r, a, s) = r) ras r] ras = 
1, A + _(r, a, s) = r) ras 7] r - as = iaa x , v4_ + (r, a, s) = 
ilrasV-ras = iraa z , and A — (r,a,s) = T] ras V-r-as = 
—iray with the standard Pauli matrices <7i (i — x,y,z). 
The quantity qp = 7m/4 is related to the deviation n of 
the average electron density from half-filling and can be 
controlled by the gate voltage. The Hamiltonian, Hq and 
'Himp = Si=i 2 'Himp, is expressed by the phase variables 
as follows, 



U »= E E|r J Ax{K- s \d x 9 ]S f+K ]S {d x cp 3& f}, 

j=p,a S=± 

(5) 

H imp=-^ / dx7 1 (x)^2raexp(~2irq F x) 

ras 

x exp{-ir((9 p+ + s9 a+ + a0 p _ + as0 a -)} , (6) 
dx^2 r ex P{— ir(2q F x + 6 p + + s9 a +)} 

rs 

x [£(x) exp {-i(0 p _ + s^ ff -)} + h.c.} , (7) 



n,2 = ZlZl 

tmp 2ira 



where K p+ = (vp+fvo)^ 1 = 1 / 1 + AV (0) / (ttv ) and 



K 



vjs/vo 



1 for the others. The Pauli matrices 



seen in Tt] mp are due to the product of Klein factors. 

With the above phase Hamiltonian, I calculate the dy- 
namical conductivity, a{u>), which is expressed by the 
memory function, M(cj), as follows Jl^ ], 



a(u) = 
M(lj) = 



-ix(0) 
u + M(oj) ' 

((F;F) LU -(F- 1 F) LU=0 )/m 

-x(o) 



(8) 
(9) 



where (A; A) u = -if dx dte^-vV {[A(x,t),A(fi,Q)]) 
with r\ — > +0, (■ ■ •) denotes the thermal average with 
respect to 7i, F = [j,TL] with j being the current opera- 
tor, and x(0) = (j; j) w =o- Since the present Hamiltonian 
conserves the total particle number, there are no non- 
linear terms including cj> p +. Then the current operator 
is expressed by j = 2v P +K p+ d x (j> p + / 'tt , which leads to 
x(0) = —4v p +Kp+/ir. To lowest order in D\ and D2, 
M(lo) is calculated as 



2nv n +K 



M(u>) 



, „ . Ki LJ 

x <^ B(—-i . 



E 



D, 



— ' (TraY 
1,2 v ' 



■ sin ■ 



TiKi (2-kT 



A", 



1 

(10) 



where K Y = (K p+ + K a+ + K p _ + K a J)/2, K 2 = (K p+ + 
K a+ + K~} +K~_)/2, cup is the high energy cut-off, and 
B(x,y) — T(x)T(y)/T(x + y) where T(x) is the gamma 
function. 

For ui — > 0, M{u>) reduces to 



M(0) - i2v p+ K p+ ~i 



Di (2-kT\ K * 2 T 2 (K l /2) 
Pi) 



(11) 
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which leads to the resistivity, 



P = a(oy 



2a ^ 

i=l,2 



/2ttT 



A", 



r 2 (^/2) 



(12) 



where 2?^ = Dia/{nVQ) = a/(ir£i) with ^ = uor^ 
being the mean free path. In the present case of 
Ki = (K p + + 3)/2, the above expression results in p = 

Pbo (2nT/uj F ) (K ^- 1 ^ 2 T 2 ((K p+ + 3)/4)/r((X p+ + 3)/2) 
where pso — 7: /2(£ 1 1 + £ 2 1 ) is the resistivity for the 
non-interacting system in the Born approximation. The 
above result shows that the contribution of the two 
kinds of backward scattering to the resistivity is ad- 
ditive. For typical nanotubes with K p+ ~ 0.2, the 
resistivity is enhanced compared to that without the 
Coulomb interaction and shows a temperature depen- 
dence as p oc T^- 1 )/ 2 ~ T-° A Q. In the pres- 
ence of the long range Coulomb interaction, the phase 
expressing the symmetric charge fluctuation, 9 p+ , be- 
comes rigid compared to the noninteracting case, which 
is expressed by the fact of K p+ less than unity. Such 
a rigid phase is easily pinned by the impurity potential. 
Thus, the long range Coulomb interaction enhances the 
resistivity. From the asymptotic behavior of Eq. ( |l0| ) for 
lu T, the optical conductivity for high frequencies is 



calculated as <t(lo) oc u> 



(K p+ -5)/2 



~ 2A , which decays 



faster than that for the noninteracting case, <t{oj) oc oj~ 2 . 
I note that the above two results of the resistivity and 
the optical conductivity do not depend on the filling, 
q F . The umklapp scattering has been known to be an 
other origin of resistivity. The umklapp scattering leads 
to p oc T ( 



2K„ 



- T-° b HI for T > v q F and 
a{uj) oc u,( 2K p+-V ~ uj~ 2 - 6 for uj > v a q F . Though the 
powers due to the umklapp scattering are similar to those 
given by the impurity scattering, it is possible to separate 
them by tuning the gate voltage. 

Next I take into account of the effects of localization by 
the renormalization (RG) method. Following Giamarchi 
and Schulz 0, averaging over the random fields, rj{x) 
and ^(x), I obtain the action, Si mp , corresponding to the 
impurity potential as follows, 



S 11 = 



X 



{■ 



2 \naj 
i(2q F x 



J Ax J drdr' 

- p+ ) cos 9 a+ sin p _ cos 9 a - 



COS' 

x cos{2qpx + 0' p+ ) cos 6' a+ sin 
+ 9 P+ ) sin 9 a+ cos 9 P 



+ sm(2qpx 

x s\a{2q F x + 9' p+ 



imp 



x sm(2qp; 

-D, (— 

\irci, 



- sin^o-- 
6»; + cos6» p _sin6i;„ 
'X + 9 p+ ) cos 9 a+ sin 9 p ^ cos 9 a . 



— 2 cos(2(7i?, 

•x + 9' p+ ) sin 9' a+ cos 9' p _ sin I 

dx / drdr'e- 1 ^"— 



-}, (13) 



sm(2qpx + 9 p+ ) cos 9 a+ cos <f) a - 



sva{2q F x + 9' p+ ) cos 9' a+ cos 
9 p+ ) sin 9, 



■ cos(2q F x 



X 



- 1 



- 0' p .,) sin 0, 



9 a + sin ( 
9' a+ sin i 



V- 



008(2(7, 

sin{2qpx + 9 p+ ) cos 9 a+ cos 

x 003(2(7, 
cos(2qpx + 9 : 



q F x + 9' p+ ) sin 9\ 



> a+ sin <j) a _ 
> p+ ) sin0 CT+ sin0 CT _ 



x s\i\(2q F x + 9' 



? p+ )cosl 



'U cos 



(14) 



1 i 0j5 = 6]S(x,t'), </>jS = <t>]8{x,T) , 

It should be noted that the intcrac- 



where 6j$ — 9jg(x,r 
and 4>' jS = 4> jS (x,T'). 
tion processes generated from the equal-time component 
in Eqs.(fl3|) and (|l4] ) are disregarded. The results given 
in the following do not qualitatively changed even when 
the such processes are included because such operators 
are less divergent than S\ mp . By calculating the various 
correlation functions to the lowest order of D\ and D 2 , 
we have the following RG equations, 



V 2 



K 



3+ 
u j+ 



K 



3- 



{3 - (K p+ 
{3 - (K p+ 



K, 



K 



p- 
-1 



K a _)/2}V U 
■K-1)/2}V 2 , 



= - {Vx/Xx + V 2 /X 2 ) K 2 + u J+ , 

= - {Vx/Xi + V 2 /X 2 ) K J+ u 2 + /l6, 

= - (V 1 K 2 _/X 1 - V 2 /X 2 ) Uj -_, 

= - (V 1 K j _/X 1 + V 2 Kr 1 /X 2 ) u 2 _/16, 



(15) 
(16) 
(17) 
(18) 
(19) 
(20) 



where ' denotes d/dl with d£ = d\n(a/a) (a is the new lat- 
tice constant), X\ = u p ^ + ' 2 u^ 

K 



-/2,.^-/2 K„_/2 y 
U„_ U„_ , A2 



p- 

and j = p or a. The ini- 
tial conditions for the above RG equations are as fol- 



v p+ /2 K 



+ /2 1/2A" P _ l/2K„ 
U p _ uj_ 



lows. 
K, 



I><(0) 



Dia/{-Kvl), K p+ (0) = u p l = K p+ , and 
ff+ (0) - K p _(0) = K a „{0) = u a+ {0) = Up _(0) = 
u a -(0) = 1. The above equations together with the ini- 
tial conditions indicate that K p _ = K a — and u p _ = u CT _. 
In addition, the RG equations are invariant under the 
transformation, T>\ «-> I? 2 ! ^ -^p— ' an< ^ -^<r- 

Therefore I can discuss the case of 2?i(0) > 2?2(0) 
without loss of generality. 

By solving the RG equations numerically and using 
Eq. (|l2|) , the resistivity is obtained as a function of the 
temperature I show the result for T>\ = 1/300 and 
T> 2 = 1/3000 for both cases with and without Coulomb 
interaction in Fig. 1. Surely, the resistivity is enhanced 
by the effects of the localization for low temperatures. 
Fig. 2 and the inset of Fig. 2 show the resistivity due 
to the inter- valley scattering, p 2 , and that by the intra- 
valley scattering, pi, as a function of the temperature. 
The quantity, pi, increases monotonically with decreas- 
ing temperature, whereas p 2 tends to zero in the low 
temperature limit. The temperature below which a de- 
crease of p 2 is observed is higher than that for the non- 
interacting case. When the inter-valley scattering is 
stronger than the intra-valley one, p 2 increases and p\ 



3 



vanishes. The result is characteristic for the localized 
regime and shows that the two mechanisms of scatter- 
ing cannot coexist. As is seen in Eqs.(6) and (^), the 
intra- valley (inter- valley) scattering pins the phases, 9 P + 1 
, and 6 a - {6 P +, &a+, 4>p-, and 4> a -)- Since the 



conjugate variables, 9 p - and 4> p -, or 6* CT _ and CT _ can- 
not be pinned at the same time, the two kinds of the 
backward scattering cannot coexist at low temperatures. 
In the present analysis, the quantitative discussion on 
the localized regime and crossover towards it have not 
done. Therefore, more detailed discussion are needed for 
understanding of disordered SWNTs with the Coulomb 
interaction. 

In conclusion, I investigated the transport properties 
of the disordered SWNTs with Coulomb interaction. I 
found that the interaction enhances the resistivity, leads 
to a power-law dependence of the resistivity as a func- 
tion of the temperature, and modifies the power of the 
frequency of the optical conductivity. In addition, it is 
shown that the intra- valley and the inter- valley backward 
scattering cannot coexist in the localized regime. Finally 
it is mentioned that the increase of the resistance with 



decreasing T observed in SWNTs |19| may be due to the 
impurity scattering studied in the present paper because 
the difference of the work function of the metallic elec- 
trode and that of the nanotube results in a downward 
shift of Fermi level in the nanotube EO) . 
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FIG. 1. The resistivity, p, normalized by pbo as a func- 
tion of 2-kT/lo f in case of Vi = 1/300 and V 2 = 1/3000. 
Here (1)((3)) is the resistivity with (without) the interaction 
derived from the RG analysis, and (2)((4)) is that with (with- 
out) the interaction derived from the perturbation theory. For 
(1) and (2), K p+ = 0.2 is used. 
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FIG. 2. The normalized resistivity due to the inter-valley 
backward scattering, p2/pB0, as a function of 2ttT/ll>f- The 
correspondence between the curves and the numbers, (l)-(4), 
and the parameters are the same as that of Fig.l. The white 
circle shows the temperature corresponding to Pi ~ 1, be- 
low which the perturbative RG analysis breaks down. Inset 
: The normalized resistivity due to the intra-valley backward 
scattering, pi/pBo, as a function of 2tvT/ujf for the same pa- 
rameters. The correspondence between the curves and the 
numbers, (l)-(4), and the parameters are the same as that of 
Fig.l. 
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